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Tensor models are a generalization of matrix models (their graphs being dual to higher- 
, dimensional triangulations) and, in their colored version, admit a 1/N expansion and a 

continuum limit. We introduce a new class of colored tensor models with a modified prop- 
agator which allows us to associate weight factors to the faces of the graphs, i.e. to the 
bones (or hinges) of the triangulation, where curvature is concentrated. They correspond to 
dynamical triangulations in three and higher dimensions with generalized amplitudes. We 
solve analytically the leading order in 1/N of the most general model in arbitrary dimen- 
sions. We then show that a particular model, corresponding to dynamical triangulations 
with a non-trivial measure factor, undergoes a third-order phase transition in the continuum 
characterized by a jump in the susceptibility exponent. 
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in ' I. INTRODUCTION 

Statistical models of fluctuating geometry are a generous source of results and ideas for physics 
and mathematics. A particularly attractive feature of many such models is that they can be 
thought as providing either a regularization or a fundamental description of quantum gravity [lj]. 
Dynamical Triangulations (DT) 0, 0j are one of the most studied examples, and have been very 
successful in two dimensions, where they are related to the large- N limit of matrix models 0, @], 
& ' and whose link to non-critical string theory in the continuum limit is well- understood 0]. 

Higher-dimensional models of DT have not been equally successful in providing a sensible con- 
tinuum limit for quantum gravity 0], leadingto degenerate geometries at large scales 0] with a 



first-order phase transition separating them [9|, llOf] . However one cannot exclude that some un- 



known essential ingredient was missing in the models analyzed so far. A non-local modification 
of DT, which goes under the name of causal dynamical triangulations or CDT [111] has produced 



substantial evidence for the emergence of an extended geometry at large scale [12j, [13|] , and signs 



of a second-order phase transition 11411 . hinting at the possibility to attain a good continuum limit. 
Recently it has also been suggested [151 ] that the effects of a non-trivial measure factor were over- 
looked in the past and could potentially lead to an improvement of the large scale behavior of DT 
models. 



Tensor models jl6l4l8l| and group field theories 19(] are the generalization of matrix models 



to higher dimensions. In particular, the colored tensor models and group field theories [2CH2 
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generate graphs dual to orientable [23[ pseudo-manifolds in any dimension. Much progress has 
been done in understanding these models. Various power counting estimates and bounds of graph 
amplitudes in tensor models and group field theories have been obtained 

0J35L 

The symmetries 

of tensor models have been analyzed either with the help of n-ary algebras 36l-l38l]. or, in a more 
quantum field theoretical approach through Ward-Takahashi identities [3^ ]. The relation between 
symmetries of group field theories and the diffeomorphism symmetry of the resulting triangulation 
has been explored 4(3, 41 ]. Solutions of the classical equations of motions 42, 43] have been derived 
and some of them interpreted as matter fields on non commutative spaces. Most importantly, for 
our purposes, the colored tensor models have been shown to possess an additional (and welcomed) 
feature as compared to non-colored models: their amplitudes are such that a 1/N expansion is 
possible, with the leading order encoding a sum over a class of colored triangulations of the D- 
sphere 2414261] . This discovery led to the possibility of new analytical investigations of DT models 



and their continuum limit in D > 3 dimensions 



27. 2 



The link between tensor models and DT leads straight to a daunting question: will any colored 
tensor model admit a richer continuum limit? The question is twofold, as at first instance one 
can wonder whether the coloring will already suffice to generate new universality classes in these 
models. The results of [27J indicate that, as far as the critical exponents are concerned, the color 
alone is not enough and the continuum limit is strongly reminiscent of branched polymerd3- The 
next question is then whether it is possible to modify the tensor models in such a way that new, 
appealing, phases would appear (like in CDT) and/or a second-order phase transition would occur 
between phases. 



The dually weighted matrix models introduced by Kazakov et al. 44| are an ideal candidate 
for introducing relevant modifications of a two-dimensional DT model via a local modification of 
the corresponding matrix model (see also 45-4^]). For example, one can impose the non-local 
condition on foliations that characterizes two-dimensional CDT just by a suitable modification of 
the propagator of a two-matrices model [481 ] . 



In the present work we extend the idea of dually weighted matrix models to colored tensor 
models (which by lack of fantasy we call dually weighted colored tensor models) and solve them 
analytically at leading order in 1/N. We then consider a particular example corresponding to a 
DT model initially proposed in 49[ and recently claimed to exhibit a new phase with promising 
geometrical properties [15(. In the large- iV limit, equivalent to weak-coupling from the gravitational 
point of view, we find a third-order phase transition and a new critical behavior. We compute 
explicitly the susceptibility exponent associated to the different phases. Although, as we are in 
a different regime^, we cannot make direct contact to the results of 15], our results show that 
a non-trivial continuum limit is possible in these models and open a small window of hope on 
non-causal DT models. 

This paper is organized as follows. In Sec. [II] we introduce the new models, discuss their DT 
interpretation and derive a Schwinger-Dyson equation relating the free energy and the connected 
two-point function of our new tensor models. In Sec. Mil we derive, in the large- N limit, a set of 
self-consistency equations which allow to exactly solve (at leading order in 1/N) the model. In 
Sec. IIVI we specialize to a particular model and study its phase transition. 



1 The reader should keep in mind however that no other characteristic of branched polymers, like say the Hausdorff 
or spectral dimension has so far been reproduced for the colored tensor models. 

2 And have not computed the spectral dimension of the emerging geometry. 
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FIG. 1. The index structure of a D = 3 vertex, in a "coarse grained" view, showing only the tensor color, 
and in a detailed view, showing the strands (indices should be read clockwise). 



II. DUALLY WEIGHTED COLORED TENSOR MODELS 

In this section we introduce a modification of the independent identically distributed (i.i.d.) 
colored tensor mode! of MM, 1244271] , which we baptize "dually weighted colored tensor model" . 




We denote Hi, for i = 0, . . . , D, the Z?-tuple of integers Hi = (nu-i, . . . , n^, riiD, ■ ■ ■ 
with riik = 1, • • • ,N. This ./V is the size of the tensors and the large N limit defined in 
represents the limit of infinite size tensors. We set riij = riji. Let ^| , ipt., with i = 0, . . . , D, be 
D + 1 couples of complex conjugated tensors with D indicefl The dually weighted colored tensor 
model in dimension D is defined by the partition function 

e N ° E = Z N (X,X) = [ d^d^e' 8 ^ 



D 



pi, Hi j 

A D i A D 

+ ^D(D-l)/i ^JlT^i + ND(D-l)/4 ^ZIT^ni ' W 
n i=0 n i=0 

^2 n denotes the sum over all indices riij from 1 to N. The index structure of a vertex is represented 
in Fig. [TJ The quadratic part is chosen such that Gaussian correlations, i.e. the propagator of the 
model, is (no summation over i, the propagator between fields of different color is zero) 

$Mft>o = - (^m> =u c ^ = u c p^ , (2) 

where C T denotes the transposed matrix. Obviously, when Cn ijPij = ^riypy the model reduces to 
the i.i.d. one. 

Graphs are made of vertices, (colored) lines and "faces". Like in matrix models, the indices of 
the tensors are associated to "strands" (the solid line in the detailed view of the vertex in figure 
[T]). A index (associated to the strand common to the half lines of color i and j of the vertex) 
is identified by the two colors, % and j. The colors are conserved along the lines, hence the "faces" 
(closed strands) of the graph are identified by couples of colors. A graph is dual to a triangulatiorJ!: 
its vertices are dual to D simplices, its lines to D — 1 simplices and its faces to D — 2 simplices. 
When computing amplitudes of graphs, one picks up the trace of the alternate product of C and 
C T along each face of a graph. We consider only matrices C which do not modify the scaling in 
1/N of the amplitude, i.e. such that limjv^oo — — is finite for all q. We do not impose for 



3 The epithet i.i.d. refers to the fact that in such models the free Gaussian measure weights independently and 
identically each component of each tensor. In the new models we are introducing this is not the case. 

4 We denote the indices of the ip field by ft. This does not denote a complex conjugation, being merely a book 
keeping device. 

5 More precisely to an abstract finite simplicial pseudo-manifold [2l| . 
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FIG. 2. The self energy E at leading order in terms of the connected two-point function G2. The labels 
0, 1, . . . D denote the colors of various lines. 



the moment any further restrictions on C, in particular C is not required to be hermitian. Using 
the counting of faces of a colored graph established in [if]] , the amplitude writes 

A(G) = (\\y( [I Tr[(C ^ )] ) N D ~T^^ , (3) 

with 2p^\ the number of vertices of the p'th face of colors ij. The number oj(G) > is called 
the degree of the graph Q [if]], and the leading order graphs are those of degree 0, which are dual 



to triangulations of -D-spheres [251 ] . Crucially for the study of the leading sector in the large N 
limit is the flowing fact [271 ]: at leading order in \/N the self energy S (i.e. one particle irreducible 
amputated two-point function) factors into the convolution of D connected two-point functions G2 , 
one for each color. Such graphs (represented schematically in Fig. [2]) are called melons^. Denoting 
g = AA, the free energy of the model is a function E(Cp n ,g) of the coupling g and the matrix 
entries Cp n . 

For C = I the i.i.d. model is in direct correspondence with DT in D dimensions, just like matrix 
models are with DT in two dimensions. Each graph is dual to a triangulation, and denoting Nk 
the number of fc-simplices, its amplitude rewrites 

where kd-2 and K£> are (with g = AA), 

kd—2 = In N, and k d = i D(D - 1) IniV - . (5) 

The grand canonical partition function of DT is given by 

Zdt = E ^) e KD - 2N °- 2 ~ KdNd , (6) 

where the sum is over all D-dimensional triangulations T, and s(T) is the order of the automorphism 
group of T. The grand canonical partition function of DT equals the free energy of the tensor model, 
Zdt = E. 

The large- N limit corresponds to kd-2 00, with kd = \D{D — 1)kd-2 ~ o m 3- ^he DT 
action can be interpreted as a Regge action for equilateral simplices, one has k>d-2 ~ 1/C thus the 
large- N limit amounts to the weak coupling limit, with vanishing (bare) Newton's constant. Note 
that this is true also for matrix models, as in two dimensions the Regge action just gives the Euler 
character of the surface, and the limit G — > is the planar limit. As in that case, in order to keep 
a finite G one would have to perform some kind of double scaling limit. 



6 This is stark contrast with the D = 2 case of usual matrix models. Indeed, the 1/N expansion in that case selects 
all planar graphs, not only the melonic ones. As shown in [2?| the results leading to the factorization of the self 
energy explicitly break down in D < 3. 
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For C/I the amplitudes associated to triangulations will be in general different, thus defining 
new models of DT. In particular, the effect of the modified propagator is to associate a weight 
factor to each (D — 2)-subsimplex of the triangulation, also called a bone (or hinge) in Regge 
calculus. As in Regge's construction it is precisely on the bones that curvature resides, the dually 
weighted colored tensor models can be a precious tool for the study of DT with more a complicated 
curvature dependence of the action. 

A. The connected two-point function 

The full connected two-point function of the dually weighted model, 

i^k) = U P ^MC) , (7) 

j 

is of course independent of the color i and factored along strands. More subtly, the contribution of 
the strand of color ij to the two-point function is in fact independent of j. This is due to the fact 
that the action is invariant under a transformation which permutes any two strands ij and ik in 
the field ip l and permutes the strands on all other fields (and also the fields between themselves) 
to restore the connectivity of the vertex. 

We are now going to show that knowledge of the full connected two-point function suffices to 
compute the derivatives of the free energy. First, we observe that from 

|E/4-(^- S H, (8) 

Hi 



it follows that 



N° ~ E W4^> + ^d- x Z = N D -UnC^P] + XN D d x E = , 

Hi, pi j j 

(9) 

which, recalling that E is a function only of g, can be rewritten as 

Xd x E = gd g E=[ [ — J ] -1. (10) 



Furthermore, we have 



'np - " . v^ np 

di,bi ■> 



which, using d a - ip (C 5a Cj) = SncC pi l + C 5 adc,- ip Cj = 0, yields 



dCnp Z^V-LJ- V °v n Pbij'JV a ^b t 

ai.bi 3 



Jpn 



(C-lpC-l^.D 



N 

Solving the model consists therefore in determining P. We will do this in the next section. 



(12) 



ti 



III. THE LEADING ORDER IN 1/N 



At leading order only melonic two-point functions contribute 27j. They are characterized, as 
we already mentioned, by the fact that the self energy of the model £ = (4>=,^p.)ipi amputated 
is given by the convolution (respecting the strand structure) of connected two-point functions. 
Supplementing this by the classical Schwinger-Dyson equation relating the full two-point function 
G2, the self energy £ and the propagator C of a field theory 

<h = 0^, (13) 

yields the system of equations 

S[(^^ft)(n*9«wj -^&h$W lpI > am ^ tat ^n c rimj)} II r " ' 

Qi 3 Ti 3 3 

amputated = gN' 2 ^ 1 £ W^A)' (M) 



Substituting the connected two-point function and performing the sums leads to 

D(D-l) m D(D-l) 

3 33 



-. T D(D-l) rn D(D-l) -,— T rr ~% — T 

I Pn m ~ 9 N — [Tr (PP T )] s- H(PP T C)n jPj = J[ Cn M , (15) 



where the index i has been erased as it plays no role. Surprising as it might be, equation (|15p can 
be solved analytically. To do so, we first introduce a matrix X = C~ 1 P, we multiply by X pj b. from 
the right and by P^ n from the left. Summing the resulting expression over all hj and pj (with 
range from 1 to N) we get 



-1— r -1— r D(D-l) m D(D-l) — r _ 

11^ = 11^ +9N~^[Tr(PP T )]^ U(P T P) ajbj ■ (16) 



3 3 3 

To solve for X, we first take aj = bj for all j and sum, obtaining 

n r-i D(D-l) rn DtD + l) 

[Tr(X)] = N D + gN-^-^[Ti(PP T )]^-^ , (17) 
and then we take aj = bj for all but one j and sum, obtaining 

r> 1 n 1 D(D-l) rr DtD+l) , m 

X a b [Tr(A)] = 5 ab N D ~ l + gN —[Tr(PP T )}^ \P T P) ab . (18) 

Combining the two equations we have 

n 1 D(D-1) r , D(D+1) , _ 

N°- 1 I + gN — 2[Tv(PP T )]^-^~ 1 P T P 
x — — > ( 19 ) 

' ^ £»(£»— 1) m D(D+1) \ n 

N D + gN-^-^[Tr(PP T )]^-^j 
and we finally get the following expression for C as a function of P, 

C = P K nm ^ i rr , C T C = V — ^ ^P T P , 20 

I + gcPW-W (l + ga^^PTpf 



7 



where a 2 = jjTv(PP T ) = -^Tr(P T P). As C T C is a function of P T P, the two commute, and 
equation (|20j) can be written as a quadratic equation for P T P in terms of C T C, whose physical 
solution is obtained by choosing the sign of the root that gives P T P = when C T C = 0: 



pTp 



1 



2g 2 a 2D(D+l)-4 C T C 



\l + ga D{D+l Y 



4ga D(D+l)-2 

\~ (1 + ga D ( D + 1 ')) 21 ^ 



-C T C 



(21) 



Finally, expanding the square root in Taylor series we can write 

q-l 



? T r=y.—( 2q ) 

2 1 ( 2 Q 
a = > -' 

^ (7 + 1 



D(D+l)-2 



l +5a D(D+l))' 



(22) 



(23) 



Equation (I22p determines P T P in terms of C T C and a, whereas (]23h implicitly defines a in terms 
of the traces of powers of C T C yielding the analytic solution at leading order in 1/N of the dually 
weighted colored tensor model. Note that the derivative of the free energy with respect to g is 





~ga D( - D+ V- 2 


q-l 

Tr[(C T C)i] 


"(i + 5a fl(^+i)) 2 V" 





gd g E 



1 

N 



D 



D(D+1) r _ j, , D{D+X) 



Tr[C _i P] -l=[l+gN —[TriPP 1 )] — 



1 



ga 



D{D+1) 



(24) 



hence studying only the self-consistency equation (j23f) for a suffices to study the critical behavior 
of the model. For C = I a straightforward computation leads from (|23p to aP = 1 + ga D ( D+1 \ 
reproducing the solution of the i.i.d. model found in 27]. 



IV. PHASE PORTRAIT OF A PARTICULAR MODEL 



In order to draw the phase portrait of a specific model it is more convenient to denote gd g E 



ga 



D(D+1) 



U and write the self-consistency equation (|23p in terms of U as 



U 



E 



ii 



2q 



(J 1 D(D+1) 



-9 



D(D+1) 



Tv[{C T C)i] 
N 



(25) 



In the remainder of this paper we deal with the model defined by a covariance C such that 

Tr[(C T C) q } = Nq~? . (26) 
Such a choice corresponds to the DT amplitude 



3 K£j-2A r D-2 — KdNd 



IL 



(27) 



where the product is over all (D — 2)-dimensional simplices (bones) of the triangulation, with qi 
being the number of D-simplices to which the bone i belongs. The DT amplitude (I27p was studied 
via numerical simulations in {49| , and more recently in 15 1 , where it was argued that a new phase 
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appears, for large enough /?, with promising geometrical properties. We will now solve analytically 
in the large- N limit the corresponding dually weighted colored tensor model defined by (|26|) . 

First note that a matrix C satisfying (|26p exists in the large- N limit. To determine C we 
can for instance diagonalize it and write N equations for the N eigenvalues, corresponding to the 
traces up to power N. These equations will always have roots in the complex domain^, thus we 
obtain in general a non-hermitian matrix C which satisfies (126|) up to q = N. Alternatively we can 
impose a weaker condition Tr \{C T C) q ] = Nq- p + 0{l) for every q > 1, and solve for the eigenvalue 
distribution in the large-iV limit, using standard techniques from matrix models. This way the 
spectrum of C T C can be chosen to be real. 

The self-consistency equation (|25p (and its physical initial condition) is now 

U = S((3, z(g, U)) , S(P, z) = f) Q 7 ) * 



9=1 



JJ 1 D(D+1) 2 

<9,U) = zj^rg 1 ^ , 9(P,0)=0- (28) 

[i + uf— 



We will denote g(/3,U) the solution of the equation (J2? 

The series S(/3, z) has radius of convergence z\, = \ for all values of /?, that is, it converges for 
all g and U under the curve 

z{g b (U), U) = - ^ g b (U) = J^^Ll . (29 ) 

* 4 2 U 2 1 

Consider first the case the case /3 = 0. We have 



U = S(0, z) = 1 2Z 2 f~^ 9(0, U) = {1 + U U)D+1 ■ 00) 

The function g(0, U) has an unique maximum in U = jj. The curves g(0, U) and gb(U) intersect 
for U = 1, irrespective of D. The reader can check that gb(U) is strictly decreasing for U < 1. The 
two curves are represented in Fig. (3h.. At the critical point U = j^,g c (0) = g(0, = ^ D _^ D+1 the 
function U becomes critical, as 

(g - g c ) ~ (U - => u - i ~ (<? - 5c )^ ^ s ~ (g - g c )^ • (31) 

This is the standard branched polymer phase, obtained in the i.i.d. colored tensor model setting 



m 



27] 



We now slowly turn on f3. The derivative of S(f3,z) w.r.t. z is = ~ ^, z )- The 



derivatives of the solution g(/3, U) of the equation (j28|) are 



P 



U = S(p,z(g,U))=>dU = dpSdp + d z S d uZ dU + d g zdg , (32) 

l-S(B-l z ) £i£±lt±^iM£i£^ 
(dg\ l-d z Sd v z_ ' D(D+i)s(/3,z)(i+s(/3,zj) 

VdUJp ~ d z sd 9 z - s^-i,^ ' 

/dg_\ dpS 
\dp)u d z Sd g z ' 



Actually it can be shown that for given N the solution is unique up to permutations of the eigenvalues. 
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S« 0.1 





(a) 



(b) 



FIG. 3. The curve g(j3, U), for D = 3 and for j3 — (a) and (3 = 0, 1, 2 (b), together with the critical curve 
z(g b (U),U) =1/4. 



where z = z(g((3, [/"),[/). The function g(/3, [/) becomes critical for = 0, i.e. at [7 = U<{(3) 

solution of 



l-c5 z 5 



g=g(P,U),U 



(33) 



P,z=z{g(fi,U),u) 

thus the critical curve g^~(/3) is characterized by 

g<((3) = g{P,U<(P)) , z<(/3) = z{g<(f3),U<(f3)) , £/<(/?) = s(p,zf{0)) . (34) 

Equation (I33p translates in parametrized form in terms of zf(f3) (dropping the argument /3 to 
simplify notations) as 



, _ q(R _ , ,< + -2- S(/3, *<)£>(£-!) 
lP ' c> D(D + l)S(P,z<){l + S(P,z<)) 

Using ([Ml) , with (|32l) and (J33]), we obtain 







dgf = fdg_\ 
dp \d0Ju 



d z SdgZ 



>0, 



(35) 



(36) 



thus the critical curve g^(f3) is increasing with (3. As (jfffz) u _ u < for all /3, the critical behavior 
remains that of (|3ip . 

However this holds only for small enough /3. Indeed we find that for /3 > /3 C > 1, g(/3, U) exits 
the analyticity domain of S(f3, z) before it can reach its first maximum. Some curves g(/3, U) for 
increasing values of (3 are represented in Fig. [3^- 

For j3 > (3 C both ^ and §jj remain finite all the way up to (and including) the boundary gb{U). 
The self consistency equation becomes 

dU\ S(0-l,z) 



D(D+l)g 



(™L\ = 

V dg )p i-SlB-1 z ) D ( D + 1 )- 2 - s ^ z ) D ( D - 1 ) ' 

' D(D+l)S(fi,z)(l+S(p,z)) 



(37) 
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g(p c .u) 



u 

FIG. 4. A schematic representation of the phase transition and the critical line g(j3, U c (/3)) (in bold), 
intersecting the line g((3 c , U). 



and the denominator stays finite. Close to z = j, the non-analytic behavior is given by S(/3— 1, z) 
(1 -4z) /3 " 1+ 5, hence 



( -3- ) ~ (9 - 9bf~ 1+ ^ => ^non-analytic ~ (5 ~ Sb)^ 4 " 

V C<7 ' non-analytic 



(38) 



Note that, as /3 C > 1, for /3 > /3 C the non- analytic part of the free energy .E is always preceded by 
an analytic part with non-zero linear and quadratic terms. 

It follows that for j3 > f3 c the critical curve g^(/3) is denned by 



g>{p) solution of z(g,S(J3,l/4)) 



1 

4 ' 



that is g^(P) = gb(S(P, 1/4)) . Along the critical curve g^{/3) we have 

U>(0) = S(/3, 1/4) , ff >(/3) = g b (U>((3)) , z>(/3) 

The derivative of (/^) computes using the implicit function theorem as 

dg> _ djjzdpS 



d/3 



dgZ 



>0, 



(39) 



(40) 



(41) 



as [/>(/?) < 1. 

The critical j3 c corresponds to the point where the two critical curves g^(/3) and g^(P) meet, 
that is 



<£(&) = <£(&), U<(P) = U>(J3), z<(P) = z>((3) 
which translates, using equation (j35|) . (|4"0]) and (j4*2|) . as the single equation 

D(£ + l)S(&,l/4)[l + S03 c ,l/4)] 



1,1/4) 



£>(£> + 1) - 2 - S{/3 C , l/A)D(D - 1) 



(42) 



(43) 



Numerically we find /3 C » 1.162 for £> = 3, /3 C » 1.216 for D = 4, and /3 C » 1.134 for D = 00. A 
schematic portrait of the transition is presented in Fig. HJ 
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The critical behavior is different at /3 = /3 C , as at g = = g c we have at the same time 



Jfe ) =0, and a non-analytic behavior of S(f3, z) with exponent /3 C + i < 2, hence 



(5 - 0c) ~ - ^o) A+ 5 U - U c ~ ( 5 - 5c ) £ ~ ( 5 - 5c ) . (44) 

In conclusion the susceptibility exponent, defined by -Enon-anaiytic ~ {g — 5c) 2-7 , is 

( \ for p < p c , 

7= < l^f for /3 = /3 c , (45) 
[ ^- 2 /3 for /3 > /3 C . 

We observe that for a large range of values of /?, i.e. for f3 < /3 C , universality holds: the 
critical exponent 7 is independent of (3. On the contrary, for (3 > (3 C we have a one-parameter 
family of different critical behaviors. A possible interpretation of such unusual behavior is that 
for j3 sufficiently large the measure term starts behaving as non-local or long-range interaction, for 
which universality is not expected to hold. 

A. Order of the phase transition 

Introducing the canonical DT partition function via 

Zdt^d, kd-2) = Yl e ~ K ° Nn 'Zdt^Nd, kd-t) , (46) 

one has that the thermodynamic limit for the DT free energy is given by 

i^oo = lim -L In Z DT)Can (N D , k d „ 2 ) ~ - hi£ c (/3) . (47) 

The order of the phase transition around (3 C is then assessed by studying the discontinuity of g c (/3) 
or its derivative at /3 C . Combining eq. (|36|) . with (JiTj) and using eq. ff33|) we obtain 

f(V = f(M. (48) 

thus the phase transition is higher than first order! 

To check the precise order of the phase transition we re-parametrize the self consistency equation 
by eliminating U 

U x ~d^+T) 2 ~ 1 f2q\ 

z = H(g,S((3,z)) , g(P,0)=0, (49) 

and denote its solution g(f3,z). We therefore have 

dz = d g Hdg + duH(dpSd/3 + d z Sdz) , 

/dg\ 1 - d z Sd v H /dg\ dyHdpS 

KdzJp d g H ' \d(3Jz d g H ' 1 ' 

The two critical curves are 

g<{P)=g(P,zf(p)) with zf{0) solution of l-d z Sd v H = , 
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gfW) =9(P, \) with z>(P) = \, i.e. g>(0) solution of ± = H(g, S(P,l/4)) 



(51) 



and meet at p c , when 



4{Pc) = - A g>(P c ) = g<(P c ). 
It follows that the derivatives of the critical couplings are 

duHdpS duH(g<(P),S{^z<(P)))d p S{P,z<^)) 



(52) 



dgf _ fdg_ 

dp \dpJz 



d g H 



dg> 
dp 



duHdpS 



d g H 



/3,l/4 



d g H\g<{p),S{p,z<(P)) 
duH(g>(p), S((3, l/4))ajg5(/3, 1/4) 
d g H(g>[fi),S{p,l/i 



(53) 



which, as we already knew, are finite and continuous at the critical point. The new parametrization 
comes in handy once we move to higher order derivatives. 

All the derivatives of g^ are finite at P = j3 c (note that H (g, U) is analytic for g > 0, U > 0, 
hence its derivatives just go to constants at /3 c ,z = |), and the derivatives of gf will differ from 



them due to terms involving the derivatives of z< 



lim 

A/3->o+ V d/3 2 



2„< 



d z g 



d 2 g> 



-A/3 



dp' 



=+A/3 



dzl 
dp 



(54) 



We thus only have to check at which order the derivatives of are non-zero and/or singular. 



In order to evaluate we use 



d(d z SduH) = 



dz< dp z Sd v H - 8 gU H§§ + d z SdpSduuH 



(55) 



dp d zz Sd v H + (d z S) 2 duuH 

The singular behavior comes from approaching the convergence radius of S(ft,z), where we have 



S(P, j — Az) « S(P, j) - 4S(P - 1, \)Az + 

4 4 4 sjir 

from which we deduce that at P = p c the first singularity is in 

dz< 



d T7 S ~ Az^ 



dp 



(56) 



(57) 



Finally we use the self-consistency equation and the expansion of S(P, z) to get Az ~ A/3 /3c 2 , 



and 



while 



dp 



A^->o+ V dp 6 



d/3 2 



(/3c) 



d*g> 



dp s 



/3 c +A/3 



2„< 



d Z Z 



dp 2 



dp 2 {Pc) ' 



2-2ff c 

A o 8? A P^° 

)- OO 



(58) 



(59) 



thus the transition is third order. 
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B. Comparison to branched polymers 

The /3-dependent behavior we found is reminiscent of that of certain models of branched poly- 



mers (BP) [50, |51[. Such models are defined by the partition function 

Z B p(ll) = Y,^ N Zbp(N) , (60) 

N 

where 

N 

z bp (n) = n^) ( 61 ) 

BP N i=l 

is the canonical partition function for rooted trees BPjv with N vertices, and n« is the degree of 
the vertex i. Making the choice 

p(n) = n~ a , (62) 

one finds a critical behavioi|f| Zbp(^) ~ (/i — /Uc) 1 " 7 , with the susceptibility exponent 7 given by 

\ for a < a c , 

f ~ 0.3237 for a = a c = 2.4787 , (63) 



Ct c — 1 

2 — a for q > a c . 

Comparing to the DT model we studied here, we see that we get qualitatively the same kind 
of behavior (up to the precise value of the critical coupling) if we identify a = (3 + 3/2. The extra 
3/2 is coming from the asymptotic scaling of the factor ^4 ~ g _3//2 4 9 , and can be understood 
as the natural entropy factor of the triangulations with trivial measure. The precise value of the 
critical coupling differs in the two cases (as for any D we have /3 c + 3/2 > a c ), but this comes as no 
surprise as critical couplings are generically non- universal quantities. For example one can adjust 
its value, without affecting anything else in (j63|) . along the lines of [EH]: introducing an additional 
weight t for the the number of vertices on the last generation of branches of the tree, one finds that 
the critical point changes with t, while the susceptibility exponent above and below the transition 
remains unaltered. For that particular modification one finds that decreasing t leads to an increase 
in a c , and that for a c > 3 the value of 7 at the critical point remains frozen at 7 = 1/2. 

It is very tempting to exploit the apparent connection to BP in order to extract other properties 
of the model. For example, a number of results are known about the Hausdorff dimension du [E^] 



and spectral dimension d$ |53j, [54| of BP. In particular, following [54| we could conjecture that, for 



7 > (that is for j3 < f3 c in our model) 

d H = ~, d s = -^—, (64) 
7 1 + 7 

while for 7 < (i.e. (3 > (3 C ) 

dji = 00 , ds = 2. (65) 

For (3 < 0c we have the standard BP phase with dn = 2, in agreement with the numerical results 
for DT [8j]. For f3 > f3 c we enter a new phase, which doesn't seem to have been observed yet in 



Note that as a result of the presence of a root, the partition function for rooted BP is interpreted as the derivative 
of a "non-rooted" model. 
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simulations of DT. From the BP point of view it is understood that for (3 > (3 C the dominating trees 
are actually short bushes, with the limit of (3 — > oo eventually represented by a bush with all the 
branches attached to a unique vertex. Given this intuition, and du = oo, one could hypothesize 
that the new phase corresponds to the crumpled phase of DT (which at (3 = is only visible at 
small K£>_2). However, as the crumpled phase has 7 = —00, such an interpretation can only make 
sense for (3 — > 00. 

Interestingly at (3 = (3 C we could get 2 < du < 00, as 7 > 0. One then could possibly tweak the 
model so that the value of (3 C at the critical point fixes either du = D or ds = 

However the reader should be warned that it is not at all obvious whether the parallel with 
BP should be extended to the Hausdorff and spectral dimensions. Sharing one critical exponent 
is of course not a sufficient reason to believe that all critical exponents are common. The relation 
between melon graphs and trees was already made evident in [23], where it was used as an exact 
bijection for the combinatorial counting. On the other hand, and most importantly, the notion 
of neighborhood in a triangulation seems to be very poorly represented by the abstract trees 
associated to melons. Nevertheless it is intriguing that a strong relation between DT (for (3 = 0) 
in the weak coupling phase and BP is supported both by numerical simulations @] (supporting 
for example du = I/7 = 2) and by theoretical arguments 55(]. Whether or not the link can be 
clarified further in the context of the colored tensor models is an open issue. 



V. FINAL DISCUSSION AND OUTLOOK 

In this paper we have introduced and solved a new class of colored tensor models. The new 
models, dubbed dually weighted colored tensor models, allow to associate arbitrary weights to 
the faces of the graphs, i.e. to the bones of the dual triangulation. Choosing the weights as in 
([26]) . we have explicitly computed the susceptibility exponent of the model, and showed that in 
the continuum limit the model admits two phases separated by a third order phase transition. 
The results are reminiscent of certain models of branched polymers, with which the susceptibility 
exponent shares the same qualitative behavior. However, we have not computed the Hausdorff 
dimension for our model, and we have at the moment no reason to believe that it agrees with that 
of the BP. 

We believe that our results are very important for the DT approach to quantum gravity. First 
of all, this is to our knowledge the first time that a phase transition in the continuum limit is 
accessed by analytical means in dimensions higher than two. Furthermore, it is the first time that 
a third-order phase transition is observed in DT, which opens up the possibility of obtaining a 



continuum limit with 2 < du < 00, without the causality condition employed in CDT [lll-ll4l]. The 
link to the numerical results of [15[ is unclear, and deserves further exploration, either by pushing 
the simulations to larger k>d-2 and larger volumes, or by trying to extend the analytical tools to 
finite kd-2- 

It should also be mentioned that the BP models we have discussed in Sec. IIVBI have been 
mapped and generalized to a balls-in-boxes models 56-58J, which has been interpreted as a mean 



field model of DT. Indeed the model successfully reproduces the first-order transition (at (3 = 
and kd-2 = ^0-2) between a crumpled and a BP phase, as in the DT simulations (9L flO I , a nd in 
the {(3, KD-2} plane it shows a phase diagram which is reminiscent of the one found in 49|] (and 



more recently revisited in 15[). In the light of that, the new phase we observed might extend at 



finite n D - 2 into a phase similar to the "condensed phase" of or the "crinkled phase" of & 



The relations (|64l) imply that branched polymers can have ds = du only for 7 = 1, thus it seems unlikely to be 
able to fit both dimensions to D simultaneously. 
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Indeed our result of a negative and /3-dependent exponent 7 seems compatible with the findings of 
0. 

In this work we have concentrated on a specific choice of weights, but the solution provided in 
Sec. IIIII can be used for any other choice of matrix C. For example, we could consider 

Tv[{C T C) q ] = Nq-^e^ 1 . (66) 

The case n > 1 is not very interesting. The modification either shifts the critical point (n = 1) or 
makes the series in ([25]) always convergent (resp. divergent) for n > 1, fx < (resp. n > 1, /i > 0). 
An interesting choice from the DT point of view is n = — 1, corresponding to the DT amplitude 

A{Q) = e KD - 2N °- 2 ~ K ° N °+^i f, 11 qf . (67) 



Such amplitude corresponds to the addition of an R 2 term to the Regge action [59[]. For n < 
(|66[) the asymptotic of the series in ()25[) are unaffected, and the large- N limit of the model would 



be the same as the one we have found for \i = 0. This is in agreement with the findings of [59f | . 
according to which the inclusion of higher-derivative terms does not significantly affect the phase 
diagram of DT. Finally, for < n < 1 we have 7 = —00 for f3 > f3 c . In such case the similarity 
between the f3 > (3 C phase and the crumpled phase of DT seems stronger than what we have 
outlined in Sec. IIVBI It would be interesting to study this model further and check whether the 
phase transition remains third-order or it becomes second or first order . 

Another option would be to modify the model by choosing C such that at the critical point 
g c the first m > 1 derivatives of the coupling become zero, rather than just the first one. In 
all likelihood this would still relate to some type of branched polymers, like the multicritical BP 



studied in [52J. 



The most important open problems in our view are to study whether the connection to branched 
polymers extends to other critical exponents, and in particular to the effective dimension of the 
triangulations. Even more important would be to access finite kd-2-, maybe via some sort of double 
scaling limit. 
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